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ABSTRACT Biomembranes are thin capacitors with the unique feature of displaying phase transitions in a physiologically 
relevant regime. We investigate the voltage and lateral pressure dependence of their capacitance close to their chain melting 
transition. Since the gel and the fluid membrane have different area and thickness, the capacitance of the two membrane phases is 
different. In the presence of external fields, charges exert forces that can influence the state of the membrane, thereby influencing 
the transition temperature. This phenomenon is called electrostriction. We show that this effect allows us to introduce a 
capacitive susceptibility that assumes a maximum in the melting transition with an associated excess charge. As a consequence, 
there exist voltage regimes where a small change in voltage can lead to a large uptake of charge and a large capacitive current. 
Furthermore, we consider electromechanical behavior such as pressure-induced changes in capacitance, and the application of 
such concepts in biology. 
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Introduction 

Biological membranes provide a barrier between cells and or- 
ganelles that serves to maintain differences in chemical and 
electrical potentials by separating molecules and ions. The 
electrical phenomena which result from transient changes in 
these electrochemical potentials provide the basis for our present 
understanding of the electrophysiology of biomembranes (T , 2). 
The cell membrane consists mainly of a lipid bilayer into which 
proteins are embedded. It is widely believed that the lipid bi- 
layer itself is impermeable to water, ions and molecules. There- 
fore, electrophysiology considers the membrane to be a capaci- 
tor, and ion channel proteins are regarded as electrical resistors. 
The nerve pulse, for instance, is considered as a propagating 
segment of charged capacitor loaded by currents through the 
channel proteins Q. 

The capacitance. Cm, defines how much charge, q, is stored 
on two capacitor plates at a fixed voltage, Vm, 



Q — Cm ■ Vm ■ 

For a parallel plate capacitor. Cm is given by 

A 



Cm — 



D 



(1) 



(2) 



with vacuum permittivity eo — 8.85410^^^ F/m and dielectric 
constant e ss 2 — 4. yl is the area of the membrane, and D is its 
thickness. 

Excitatory processes in cells are typically accompanied by 
changes in voltage. During nerve pulses, for instance, the volt- 
age changes transiently by about 100 mV in one millisecond. 
Voltages as large as 100 mV are also typical in the voltage 
clamp experiments that are used to measure protein conduc- 
tances (4). In electroporation experiments, used to transport 



drugs or DNA into cells, voltages can be on the order of several 
100 mV (i5i|6). Stimulation voltages for nerve pulses can be of 
order up to 1 V (e.g., (|7]l). 

A change in voltage leads to a capacitive current because 
the charge on the capacitor changes. The capacitive current 
induced by a change in voltage is given as 
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In electrophysiological models such as the Hodgkin-Huxley 
model for nerve pulse propagation and in the interpretation of 
voltage clamp experiments, it is assumed that the capacitance of 
biomembranes (in particular of nerves) is independent of volt- 
age (i.e.. Cm is constant) so that the second term on the right 
of this equation is zero (e.g., (l2j|3jl). This is equivalent to as- 
suming that membrane dimensions are unaffected by electri- 
cal phenomena and that the excitation of membranes does not 
change their dimensions. The second of these assumptions is 
known to be incorrect because changes in the thickness of nerve 
membranes during the action potential have been observed (8j- 
fT4l. Furthermore, there are numerous reports in the literature 
on voltage induced changes in membrane bending, i.e., caused 
by flexoelectricity or mechanoelectricity JTSl - frTt . In the re- 
cent years, we have proposed that the voltage changes during 
the nerve pulse are actually related to changes in capacitance 
a3,J^-i22j|. 

In this theoretical paper we show that the assumption of 
constant capacitance is incorrect, especially if one is close to 
chain melting transitions in the lipid membrane. Biological 
membranes display transitions close to physiological temper- 
ature. Heat capacity maxima are typically found 10-15 degrees 
below physiological or growth temperature, both for bacterial 
membranes from E.coli and bacillus subtilis, for lung surfac- 
tant (13) and for nerves from the spine of rats (S. B. Madsen, 
N. V. Olsen, unpublished data). The fluid state membrane is 
thinner than the gel membrane. Simultaneously, the area of the 
fluid membrane is larger (23). This implies that the capaci- 
tance of the fluid membrane is larger than that of the gel mem- 
brane. Thus, any phenomenon in the biomembrane related to 
a phase transition will influence its capacitance. For instance, 
the temperature of the phase transition is influenced by charges 
on the capacitor because electrostatic forces act on the capac- 
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itor plates. Further, hydrostatic pressure (l24b . lateral pressure 
(fT3l l or the addition of drugs like anesthetics (|25]) influence the 
position of the melting transition. Therefore, the capacitance 
must be also a function of voltage, hydrostatic pressure, lateral 
pressure, and the concentration of anesthetics. 

We show here that close to a transition voltage is able to 
change the dimensions of a membrane and the capacitance. We 
also show how the electrical properties of the membrane are 
affected by the application of lateral pressure or tension in a 
membrane. This gives rise to pressure-induced capacitive cur- 
rents or pressure-induced voltage changes. While we restrict 
ourselves primarily to the phenomenon of electrostriction, i.e., 
the force that capacitive charges exert on the capacitor, we also 
discuss polarization effects. 



Theory 

Consider a capacitor whose equilibrium properties depend on 
voltage Vm only, i.e., all other intensive variables of the system 
such as pressure and temperature are kept constant. We write 



\dVrr, 



(4) 



Here we introduce the function C„i = (dq/dVm), which we 
call the capacitive susceptibility. Note, that the definition of 
Cm differs from that of the capacitance that is given by Cm ~ 
q/Vm- According to eq. eq. Q corresponds to 
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This equation takes into account that the changes of the charge 
on a capacitor are not solely due to voltage changes but also to 
voltage-induced changes in capacitance, which are described 
by the capacitive susceptibility 



C — C 



Vm 



dC„ 
dVm 
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capacitor plates. In the absence of other forces, increasing the 
voltage on a capacitor will generally tend to deform the capac- 
itor such that its thickness is reduced. 

Consider a membrane with fixed thickness D, a capacitance 
Cm (given by eq. (j2|) and a transmembrane voltage, Vm. The 
field across the membrane is E — Vm/D (assuming a uniform 
dielectric constant in the membrane interior), and the charge is 
q = Cm. Vm. The force, acting on this capacitor is given by 

(Ha 

^-'-E-q^'-'^q^l^. (7) 
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The force T and the field E are vectors normal to the mem- 
brane surface. Assuming that the fluid membrane has a capac- 
itance of Cm ~ 0.5 ^F/cm^ = 0.5 • 10^^ F/m^, this results 
in a pressure on the membrane of p ~ lO^N/m^ at 100 mV, 
which corresponds to 0. 1 bar. Due to its quadratic dependence 
on voltage, this pressure is 100 times larger {p = 10 bar) at 
Vm = 1 volt. In contrast to hydrostatic pressure this pressure 
has a direction normal to the membrane. This implies that in- 
creasing this force results in a reduction of thickness and an 
increase in area. Since the melting of a membrane is linked to 
an increase in area, an increase of transmembrane voltage can 
therefore potentially melt a membrane. 

At constant voltage, the work done by the electrical field 
upon melting of the membrane is given by 



1 A 
TdD = -eoV^ / e — dD , 

J Dn 



(8) 



where Dg and Df are the thickness of the gel and the fluid 
membrane, respectively. For constant area this relation yields 
the familiar expression for the work done on a capacitor, ^CmV^- 
However, the membrane area does not stay constant here. 

The area in the fluid state of DPPC is 24.6% larger than 
in the gel state and the thickness is 16.3% smaller (23). If we 
assume a dielectric constant e independent of voltage and mem- 
brane state, this difference is 



^ , Ag. (1 + 0.246) 
" \Dg • (1 - 0.163) 



= 1.49 cr 



(9) 



If the capacitance Cm is independent of voltage, we obtain 
Cm = Cm - However, the last term on the right hand side of 
eq. (j6]l can become large close to transitions in biomembranes 
as we will show below. In the context of transitions, this term 
can be considered an excess capacitance. It is proportional to 
the voltage. Therefore it is zero at zero voltage where the ca- 
pacitor is not charged. In contrast, the function C,„ has a finite 
value since it depends only on the dimensions of the capacitor 
The capacitive susceptibility Cm is fully analogous to other 
susceptibilities such as heat capacity, {dH/dT)p, isothermal 
volume compressibility, —{dV/dp)T, and isothermal area com- 
pressibility, — {dA/dlV)T. The capacitive susceptibility has been 
used before, e.g., by {26h 

Electrostriction 

Electrostriction is the generation of a mechanical force on a ca- 
pacitor by the electrostatic attraction of the charges on the two 



Thus, the capacitance in the fluid phase is about 1 .5 times larger 
than that of the gel phase, and the force across the membrane 
is 1.78 time larger than in the gel phase. In the presence of a 
voltage difference, a sudden change in membrane state can lead 
to significant capacitive currents (see discussion). 

The assumption made above of a voltage- and temperature- 
independent dielectric constant may not be correct. Paraffin oil 
has a dielectric constant of 2.2 — 4.7 and olive oil has 3.1 while 
paraffin wax has 2.1 — 2.5 (28). Therefore, it may be possible 
that the dielectric constant in the fluid membrane is somewhat 
higher Lacking reliable data for membranes, we consider e to 
be constant. 

Voltage dependence of the melting temperature 

In the following we describe the melting of a membrane in 
the presence of a transmembrane voltage. We approximate 
the excess heat capacity in the absence of voltage, Acp.o — 
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d{/S.Ha{T)) / dT, by assuming that the transition is governed 
by a van't Hoff law, i.e., by a two state-transition from gel to 
fluid with a temperature dependent equilibrium constant K{T). 
The temperature dependence of the excess enthalpy, AHq{T) 
is given by 



AiJo(T) 



K{T) = exp 



KiT) 



K{T) 
AHo 
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(10) 



where n is a cooperative unit size that determines how many 
lipids undergo a transition at the same time. For DPPC, the 
total excess enthalpy of the transition is AHq = 39 kJ/mol, 
and the melting temperature is r„i — 314.2 K. The transition 
of DPPC unilamellar vesicles is reasonably well described by 
n = 100. (For details of this calculation see (l29T l.) We de- 
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Fi gure 1 1 The influence of electrostriction on the melting of mem- 
branes. Top: Excess heat capacity of DPPC LUV at three different 
voltages. Bottom: Change of the transition temperature of DPPC 
as a function of voltage (see eq. \15\). 



note the membrane volume as ^(r) — Vg + AV{T) and the 
membrane area as A{T) ~ Ag + AA{T), where Vg and Ag are 
the volume and the area of the gel state lipids, and AV{T) and 
AA{T) are the temperature dependent changes due to melt- 
ing. In previous publications we have shown that the changes 



in both volume and area are proportional to the excess enthalpy 
during the chain melting transition (.23ll24l l: 



AV{T) ^ j^AHoiT) 
AA{T)^j^AHo{T) , 



(11) 



where 7^ = 7.8 • lO^^^m'^/J is a constant that is practically 
independent of the lipid species or the lipid mixture C24j , and 
■y^ = O.SOm'^/J. We assume that a similar relation holds for the 
mean thickness. 



(12) 



with 7^ — —2.49 • lO^^^m/J. The thickness of the gel phase 
membrane, Dg, is 4.79 nm (for DPPC), and it is Z?/ = 3.92 nm 
in the fluid phase. DPPC in the gel phase has an area of Ag = 
0.474 nm^ per lipid and a membrane area 1.43 • 10^ m^ per mol 
of lipid (23). Here and below we will assume that the above 
proportionality to the enthalpy is valid. Further, we assume 
that the voltage dependence of the pure lipid phases is small. 
According to eq. (Is]), the enthalpy change, AH{T), at constant 
voltage, V,n, of the membrane at temperature, T, is given by 



Ai^(K^,T)=Ai^o(^) + ieo v;2/^ 



(Ag+7^Ago(T)) 



(Dg+7DAffo(T)) 



'^-fodAHoiT) , 



(13) 



Making use of (1 + x) 
constant e, we finally obtain 

Ai7(K„,r) = Ai/o(T)- 



1 — 2x for small x and assuming 



(14) 
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where AHo{T) is the temperature-dependent enthalpy in the 
absence of voltage described by eq. ( [T0| . One can now deter- 
mine the temperature dependence ofA^H{Vm,T). For T 3> 
Tm, AH{Vm,T) assumes a constant value, which is the ex- 
cess heat of melting. It is a quadratic function of voltage: 
AH{Vm) = AHq + aoVl with uq = -141.7 [J/V^] using 
e = 4. 

AHq is the melting enthalpy of the membrane in the ab- 
sence of voltage with melting temperature Tm.o = AHq/ASq, 
and AH{Vm) is the melting enthalpy of the membrane with 
voltage Vm. The transition temperature T„i in the presence of 
voltage can be written as 



T 



AH{Vm) AHo + aoV,l 



ASo 

= Tm.o + 



ASo 



ASo 



aVji)T„ 



(15) 



with a = ao/AHo = -0.003634 [l/V^], and ASo = 
AHo/T„ifi- The melting profiles and the melting tempera- 
ture as a function of voltage are shown in Fig{T] One obtains a 
shift of the transition temperature towards lower temperatures 
of -11.4mK for Vm = 100 mV, of -1.14K for Vm = 1 V and 
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of —114 K for V„i =10 V. The effect of electrostriction on the 
melting temperature is obviously small at physiological volt- 
ages. However, it is large for voltages of more than 1 V. 
The above calculation applies to a symmetric membrane. Ac- 
cording to the results of ( 30 ) the quadratic voltage dependence 
of T,n can be shifted on the voltage axes when the membrane 
is asymmetrically charged. In this case, the maximum melting 
temperature can be achieved for voltages different than zero. 

Voltage and temperature dependence of the ca- 
pacitive susceptibility 

It has been shown that outside of transitions, capacitance changes 
induced by voltage are small (31-33) but quadratic in voltage. 
As above we will therefore assume that the voltage dependence 
of the capacitance in the transition regime is much higher than 
that of the pure phase. I.e., we assume that Ag and Af as well 
as Dg and Df display a negligible voltage dependence. The 



given by 



f,r ^^ Ag,i + AAiVra,T) 

'^^-^^^-''"Dg,, + ADiV„,,T) 



(16) 



and at voltage Vm + dV by 



_^ ^^ Ag,i+AAiV^ + dV^,T) 
Dgei + AD(Vm + dVm, T) 
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Figure 2: Bottom: Area and thickness changes as a function of 
voltage at T = 311 A". The temperature, T — 311 A", is below the 
melting point of DPPC in the absence of voltage. Top: Capaci- 
tive susceptibility, Cm, of DPPC LUV as a function of voltage. The 
shaded area indicates the excess charge of the voltage-induced tran- 
sition. The dashed line is the voltage dependent capacitance. Cm,. 



charge on a membrane at temperature T and voltage Vm is then 



(K^+C^Kn) 

(17) 

where and AD are again proportional to AH (cf., eq. 
([T4]i). For fixed voltage, the capacitance is a function of tem- 
perature only. For fixed temperature, it is a function of voltage 
only. The area and thickness of the membrane at T = 311 K are 
shown as a function of voltage in Fig. |2](left). The capacitive 
susceptibility. Cm, is now given by 



dq q{Vm + dVm.T)-q{Vm,T) 



dV,r. 



dVr. 



(18) 



If the temperature T is below the melting temperature, Tm,vi, 
voltage can induce a transition in the membrane. Cm displays 
a maximum at a transition voltage that depends on the experi- 
mental temperature. This is shown in Fig. |2](right) for the case 
of T = 31 IK. The units are given in absolute units (per mol of 
lipid) since the area of the membrane is not constant. For com- 
parison, the specific capacitance in the gel state is 0.74 jjFlcvc?, 
and in the fluid state it is 0.93 /iF/cm^. 
The charge on the capacitor is given by 



q{Vm.T) 



(19) 



The charge as a function of voltage is shown in Fig. |3](left) for 
various temperatures. One can see that the charge undergoes a 
stepwise change at the transition voltage. We call this change 
in charge the "excess charge". Ago- It is given by (cf., eq. (|6]l) 



Vm 



dVm 



dV„. 



(20) 



where the integral is from a voltage below to a voltage above 
the transition. The excess charge corresponds to the shaded 
peak area in Fig. |2](right). It has a value of 634 C/mol. 

One can also determine the capacitive susceptibility. Cm, 
and the capacitance. Cm, as functions of temperature. This is 
shown in Fig. |3] (right) for several voltages. For Vm = 0, the 
capacitive susceptibihty and the capacitance are identical. 

Fluctuations 

The heat capacity of a membrane at constant pressure is the 
temperature derivative of the mean enthalpy and is given by 



d{H) 
dT 



with 



(H) = 



E,g»cxp(-g,/fcr) 

E,exp(-iJ,/fcT) 



(21) 



where (H) is the statistical mean of the enthalpy averaged over 
all possible microstates of the system with enthalpy Hi. The 
enthalpies of the microstates are given by: 



H, = E,+ pVi + HA, + ^qi 



(22) 
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Figure 3: Left: Voltage-induced transition of DPPC LUV at four different (constant) temperatures. Top panel: the charge as a 
functions of voltage. Bottom panel: capacitive susceptibility as a function of voltage. Right: Capacitive susceptibility Cm as a 
function of temperature for five different voltages. The dashed lines represent the capacitance Cm as a function of temperature. 



with the intensive variables pressure p, lateral pressure H, elec- 
trostatic potential 5* and the conjugated extensive quantities, 
internal energy Ei, volume Vi, area Ai, and charge qi. 
Eq. (21 1 immediately leads to 



d{H) 
dT 



; p,n, ... = const. (23) 



which is one of the fluctuation relations. Similar relations are 
found for the specific isothermal volume compressibility k^, 
and the area compressibility (I23l l: 



j_(d{v)\_{v')-{vy 

{V) { dp 
1 fd{A) 



similar susceptibility and can be written as 



dq 



dq _{q')-{qr 



dV„, 



kT 



T,p,Il = const. 



(26) 

All of the above susceptibilities are derivatives of exten- 
sive variables with respect to the conjugated intensive variables. 
For such susceptibilities, the fluctuations (^^) — are 
quadratic forms and therefore always positive. Heat capacity, 
volume and area compressibility, and capacitance must always 
be positive definite functions. For this reason, one also finds 
that the integrals of the susceptibilities are always be positive: 



(A) \ dU 



(V) kT 

{A') - {AY 
(A) kT 



; r,n,* = const. (24) AH 



T,p,-^ = const. (25) 



Like the heat capacity, these compressibilities have maxima in 
the melting transition (1231) . 

The capacitive susceptibility of the membrane. Cm, is a 
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for 
for 
for 



T2 >Ti 

P2 > Pi 

n2>ni (27) 

Vrn,2 > Kn.l ■ 



This implies that an increase in voltage across a membrane 
must result in an increase in charge so long as intensive vari- 
ables other than Vm are kept constant. This is in agreement 
with the findings in Fig. [3](left). 
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For derivatives of extensive quantities with respect to non- 
conjugated quantities, the fluctuations are no longer positive 
definite forms, and negative values can be obtained. For in- 
stance, the volume expansion coefficient is given by d {V) jdT 
=[{VH) - {V) (H)] /kT^. For water at 0°C it is negative. 
Similarly, derivatives of other extensive variables with respect 
to non-conjugated intensive variables, such as d (q) / dT or 
d (q) /dH, can also be negative (see below). 

Piezoelectricity 

A material is said to be "piezoelectric" if the application of a 
force produces an electric field (and vice vera). "Piezo" orig- 
inates from the Greek word for pressure, and we will there- 
fore use the term 'piezoelectric' as synonymous with 'elec- 
tromechanical' in the sense of pressure-induced voltages across 
membranes. At fixed temperature. 
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(28) 



where T is the force normal to the membrane. Since thickness 
changes in the melting transition are coupled to area changes, 
this leads to the relation 
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where n is the lateral pressure of the membrane. In the fol- 
lowing we will focus on lateral pressure changes. At constant 
voltage, the enthalpy change, AH{T), of the membrane at tem- 
perature T due to a lateral pressure is given by a modified ver- 
sion of eq. ( [T3] l: 

AH{Vra,T,U) = AHo{T) + AWciVrn) + AWAili) 

AWa(JI) = IlAA = Ilj'^AHa{T) , (30) 

where AMOi(n) is the work done to change the area of the 
membrane from Ag to A. With the help of eq. ( 14 1 this leads to 
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and the melting temperature is given by 
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" " ASo ~ ^ 
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-0.003634 [W^] (see eq. (15i, and 7 



(31) 

(32) 
0.89 



with a 
m^/J. 

The charge on a membrane at temperature T, voltage 
and lateral pressure n can be written as 



q{Vm,T,Il) = eeo 



Ag,i + AAiVrn,T, n) 
Dg,i + AD{V^,T, n] 



(33) 
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Fi gure 4: Changes in electrical properties induced by lateral pres- 
sure changes. Top: Change in the charge on a capacitor upon 
changes in lateral pressure at fixed voltage ofVm =1 V- Bottom: 
Change in voltage induced by lateral pressure changes at fixed 
charge on the capacitor (Cm= 1430.7 C/mol corresponding to the 
charge on the fluid phase membrane at Vm=l V). Increasing lateral 
pressure renders the membrane more solid. 



where A A and AD are again proportional to AH (cf. eq. (31 1) 
Now we can determine (1) how the charge on the capacitor 
changes with changes in lateral pressure at constant voltage and 
temperature (and vice versa) and (2) how the voltage changes 
at constant temperature and constant charge with changes in 
lateral pressure (and vice versa). 

Fig. 4](left) shows case 1 for a fixed voltage of Vm= 1 V and 
T= 315 K. Under these conditions, the membrane is in the fluid 
state when the lateral pressure n is zero. Increasing pressure 
renders the membrane more solid and the capacitance of the 
membrane decreases. This leads to a release of charge from the 
capacitor (i.e., pressure induced capacitive currents). One can 
define the corresponding susceptibility 



[iv = 



dn 



(34) 



which is shown in Fig. |4](right). This susceptibility represents 
the change in the charge on the membrane due to an increment 
in the pressure at constant voltage. Since it is a derivative of the 
extensive variable q with respect to the non-conjugated exten- 
sive variable n, it can have a negative value. 
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Figure 5; Top: Change in voltage induced by lateral pressure 
changes at fixed charge on the capacitor (Cm = 1430.7 C/mol cor- 
responding to the charge on the fluid phase membrane at Vm=l V)- 
Increasing lateral pressure renders the membrane more solid. Bot- 
tom: The corresponding susceptibility, Pq = {dVm/dH) . 



When the charge in eq. ( |28] l is kept constant {dq — 0), we 
obtain 

(In), 



dn = Pqdn . 



(35) 



This is shown in Fig. [5] (left) where the charge is fixed at 
q =1645.4 C/mol, which is the charge on the fluid membrane 
at a voltage of Vm = 1 V. The work done on the membrane by 
pressure is converted into work done on the charges by chang- 
ing their distance. Increasing pressure at fixed charge leads to a 
larger voltage across the membrane. We call this a piezoelectric 
effect. The coupling constant. 



fdV„, 



\ an 



(36) 



is shown in Fig. [5] (right). It is identical to the term in rectan- 
gular brackets in eq. p5| ). A situation of constant charge could 
be present in a membrane during a sudden reversible change in 
lateral pressure when charges have no time to dissociate. 



Discussion 

Lipid membranes are thin capacitors that are unique in their 
property to display transitions in capacitance. Here, we have 
offered a theoretical framework for describing the capacitance 
and the capacitive susceptibility of membranes in the transi- 
tion regime. Since biomembranes are close to such transitions, 
this is of immediate biological relevance. Charges on mem- 
branes create forces on the membrane that tend to compress 
the membrane normal to the membrane surface. This effect is 
called electrostriction (27). It leads to a lowering of the melting 
temperature of the membrane. Therefore, a change of trans- 
membrane voltage can induces melting transitions. As a conse- 
quence, there is an excess charge linked to the voltage-induced 
transition. This is expressed by the capacitive susceptibility. 

Cm = dq/dVm, which displays a pronounced maximum at the 
transition. This behaviour is very similar to that of the heat 
capacity as a function of temperature and the volume and area 
compressibilities as a function of hydrostatic and lateral pres- 
sure. These functions are all derivatives of extensive thermody- 
namic variables with respect to their conjugated intensive vari- 
able. Using the fluctuation theorem, we have shown that the 
capacitive susceptibility is related to fluctuations in charge and 
must therefore always be positive. The maximum of the capac- 
itive susceptibility describes the fact that in particular voltage 
regimes close to transitions a small change in voltage may lead 
to a large uptake of charge, i.e., the membrane is very suscepti- 
ble to small variations in conditions. 

We further show that there exist electromechanical cou- 
plings in lipid membranes. One finds a straight-forward con- 
nection between lateral pressure and the charge on a membrane 
(at constant voltage) that also assumes a maximum in the melt- 
ing transition of the membrane. This implies that changes in 
pressure can give rise to capacitive currents. This effect is de- 
scribed by the derivative of an extensive variable (charge) with 
respect to a non-conjugated variable (lateral pressure). Such 
susceptibilities may have either positive and negative values. 
At constant charge, there also exists a coupling between volt- 
age and lateral pressure. This effect corresponds to piezoelec- 
tricity. The coupling constant assumes a maximum in the melt- 
ing transition. This effect is very important for excitatory pro- 
cesses (for instance nerve pulse) and is discussed below. In the 
present paper, we have not considered other potentially signifi- 
cant effects such as field-induced changes in polarization. Such 
effects are relevant and clearly exist in lipid monolayers. This 
will be the focus of a future publication. 

Due to their different capacitance, gel and fluid membranes 
carry different charges when a constant voltage is applied. An 
interesting consequence is that the presence of a transmem- 
brane voltage will lead to charge separation in the melting regi- 
me where gel and fluid domains coexists. If the membranes 
contain a fraction of charged lipids, this will lead to an asym- 
metric accumulation of charged lipids in the fluid domains. 
Thus, the shape of phase diagrams will be influenced by the 
application of an electrical field across the membrane. Not sur- 
prisingly, besides temperature, pressure and the concentrations 
of the components, the transmembrane voltage is also a vari- 
able that determines phase diagrams and should be included in 
Gibbs' phase rule. 

The influence of voltage on the capacitance of black lipid 
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membranes had been investigated both theoretically and ex- 
perimentally in several studies in the 1970's (l30l [311 134l - l36t . 
All of these studies assume that electrostriction is the domi- 
nant effect leading to a reduction in membrane thickness and 
an increase in area due to an increase in voltage. Those studies 
were done far away from the phase transitions temperature of 
black lipid membranes. In agreement with our derivations here, 
they found that the dependence of the capacitance should be a 
quadratic function of voltage, i.e.. Cm oc (1 + aV^), where a 
is a constant. This constant is strongly influenced by the pres- 
ence of solvent in the black lipid membranes (32). Solvent-free 
membranes display a much smaller voltage dependence of the 
capacitance on voltage. In an interesting experimental study, 
Alvarez and Latorre ( 30 ) found that changes in capacitance in 
asymmetric membranes are influenced by a resting potential, 
Eq, such that Cm oc {1 + a{V — Eq)'^). Alvarez and Latorre 
discussed this finding in the context of nerve pulse propagation 
and the measurement of gating currents. They suggested that 
the membrane itself could display capacitive currents similar to 
gating currents. The novel aspect of the present paper is the 
recognition of the profound effect that the melting transition 
has on the nonlinear behavior of the membrane capacitance. 
While previous authors have assumed a constant compressibil- 
ity of the membrane, we have made use of the fact that the 
compressibility is dramatically increased in transitions. 

The two terms of the capacitive current Cm ■ dVm/dt and 
Vm ■ dCm/dt in eq. ^ can have different time dependences. 
The first term is fast and is largely determined by the electri- 
cal resistance of the aqueous medium. The second term has a 
time scale given by the relaxation time of the membrane, which 
is slow in transitions. It should therefore be possible to distin- 
guish the two contributions to the capacitive current experimen- 
tally. We have previously shown that relaxation time scales can 
be as large as 30 seconds for artificial membranes at the tran- 
sition maximum, and we estimate time scales of up to 100ms 
for biomembranes (37, 38 ). Besch and collaborators investi- 
gated HEK293 cells and found currents with timescales in the 
10ms regime after voltage change (39 ), which indicates that 
these currents are related to time-dependent changes in mem- 
brane geometry rather than being caused by charging a mem- 
brane with a constant capacitance. Slow capacitive currents 
after voltage changes were also found in rat nerves (40) but in- 
terpreted as gating currents. 

The electromechanical aspect of membranes has received 
considerable interest in the past. Ochs and Burton showed in 
1974 that black lipid membranes made of egg-PC and choles- 
terol show electromechanical behavior under voltage clamp con- 
ditions ( flSl l. They applied an oscillating pressure difference 
across the membrane and recorded capacitive currents that they 
described by Ic = Vm ■ dCm/dt (cf., eq. (j3]l). 

As mentioned above, the interesting issue of polarization 
has not been treated here. However, much of the literature 
on electromechanical phenomena in membranes is dedicated 
to polarization caused by membrane curvature. R. B. Meyer 
proposed in 1969 that liquid crystals of molecules that possess 
electrical dipole moments should also be piezoelectric (AV). In 
particular he stated that "another possible application of these 
effects may be in interpreting some of the potentials and ion 
distributions observed in liquid-crystalline biological structures." 
Following this idea, Petrov and collaborators proposed that mem- 



branes should also display piezoelectric behavior (^T). Lipid 
monolayers have significant dipole moments due to the polar 
head groups and the oriented associated water layer. For a sym- 
metric membrane of a zwitterionic lipid, the polarizations of 
the two layers should cancel because they have opposite ori- 
entations. However, in the presence of curvature, this sym- 
metry argument does not apply, and one expects an electrical 
field generated by curved membranes. Petrov named this phe- 
nomenon "flexoelectricity" and demonstrated the effect in ex- 
periments similar to those by Ochs and Burton (15) but inter- 
preted differently as curvature-induced polarization (l43l l44l) . 
In several further papers the authors applied the concept of 
flexoelectricity to biomembranes and proposed a coupling of 
flexoelectricity to ionic currents through channel proteins 6451 
46 ). Since lipid membranes in the absence of proteins also dis- 
play ion-channel-like characteristics close to transitions (l47l - 
i49j), Petrov's considerations are also valid here. It is reason- 
able to expect that the polarization, P, of membranes changes 
significantly in the phase transition regime and that the related 
susceptibility, dP/ dE (with E being the electrical field), should 
have a maximum in the transition. This must be so since the 
dipole moment of a gel monolayer differs from that of a fluid 
monolayer. Upon membrane bending one should find an asym- 
metric distribution of gel and fluid lipids in the two oppos- 
ing monolayers (50) leading to an effective polarization of the 
membrane. Unfortunately, we are not aware of any reliable ex- 
perimental data or theoretical estimates of the magnitude of this 
polarization. To our knowledge, the above authors did not in- 
vestigate the effect of the phase transition on flexoelectricity. 
However, polarizations induced by curvature are completely 
analogous to voltage changes induced by lateral compression 
as described above by eq. (35 1 and Fig. |5] Interestingly, 
Helfrich investigated the effect of voltage on the phase tran- 
sition temperature of 3-dimensional liquid crystals as early as 
1970 (BTl) . He found that the shift of the transition is given by 
Ar„j = ^Tmfi£oAeE^/AHQp, where Ae is the difference of 
the dielectric constant between liquid and solid phase, p is is 
the density and E is the electric field. Assuming E/D = V„, 
this law is analogous to eq. 
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It is known experimentally that Upid monolayers have large 
dipole potentials of order 300mV — somewhat higher for gel 
than for fluid phase monofilms (32115311. This is frequently at- 
tributed to the dipolar nature of the lipid head group and the 
associated water One can influence the state of lipid monolay- 
ers in experiments by an applied field. At positive voltages a 
liquid-expanded monolayer becomes more solid, while the op- 
posite effect is observed when the field is reversed. This rules 
out the possibility that the influence of voltage is due to elec- 
trostriction (PhD Kasper Feld, NBI 2012), which is indepen- 
dent of the direction of the field. For membranes, this is less 
clear. Antonov and collaborators (54) measured the voltage 
dependence of the lipid chain melting transition via the effect 
of voltage on membrane permeability, which displays a maxi- 
mum at the melting point. For synthetic black lipid membranes 
made of either DPPC or DPPA (dipalmitoyl phosphatidic acid), 
they found an increase in melting temperature that was well de- 
scribed by 



r„(DPPC) 

T„,(DPPA) 



= 2,lhA[K] 
= 332.8[is:] 



2<d.b[K/V] ■ Vm 

bb.l[K/V]-Vm 



01) 
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This corresponds to a linear shift of T„i of +1.03 K at Vm = a velocity close to the speed of sound with many similarities to 
50mV for DPPC and a shift of 2.75 K for DPPA towards higher the nervous impulse. In fact, Tasaki and collaborators showed 
temperatures, respectively. This effect is opposed to the trend in various publications that such mechanical pulses exist (|8]- 
towards lower temperatures predicted above based on electrostric- fT2ll60b . Further support comes from atomic force experiments, 
tion, which is in agreement with previous predictions from Sugar that show mechanical signals in synapse in phase and propor- 
tional to voltage changes ( 14). Further, light scattering tech- 
niques demonstrate that voltage changes are accompanied by 
changes in nerve dimensions ( |6T1I62] |. 

One might reasonably assume that polarization pulses can 
also propagate in membranes. In lipid monolayers, such pres- 
sure pulses have recently been demonstrated experimentally by 
Schneider and collaborators (63). The pressure pulses are ac- 
companied by voltage pulses that are directly related to and in 
phase with the pressure pulse (M. F. Schneider, private com- 
munication). The coupling displays a maximum in the phase 
transition regime in an agreement with the concepts proposed 
here. Such experiments are important and will eventually lead 
to a complete thermodynamic picture of the capacitive suscep- 
tibility and the electromechanical behavior of biomembranes 
and nerves. 



However, Antonov's finding is remarkably close to the 
calculation of Cotterill (56) made by considering the polariza- 
tion of the monolayers. The shift of Tm by Antonov was found 
to be linear within experimental accuracy as expected from Cot- 
terill's calculation. We do not believe, however, that the deriva- 
tion of Cotterill is theoretically sound. 

An increase of the transition temperature with increasing volt- 
age poses an interesting theoretical problem. Under such cir- 
cumstances, the fluid membrane can be made solid by voltage 
leading to a lower capacitance. This would render the excess 
capacitive susceptibility in the transition negative an would re- 
sult in a negative excess charge. According to the considera- 
tions of fluctuations above, this can hold only if there are cou- 
plings between the charge on the membrane and non-conjugated 
intensive variables. It remains to be seen whether the simulta- 
neous change of capacitance and polarization allows for such 
behavior However, if this were the case, an increase in voltage 
could possibly result in a decrease of charge, i.e., in capacitive 
currents against the applied field. It should also be noted that 
the publication of Antonov and collaborators (54) is the only 
one on voltage-induced shifts in transition temperature that we 
are aware of, and independent experiments verification would 
be useful. In this respect, the shift of the quadratic dependence 
on voltage due to a pre-existing polarization found by Alvarez 
and Latorre (30) might contain the answer to the present prob- 
lem. This suggests that the membranes in the experiments of 
Antonov were not symmetric. Additional studies of polarized 
membranes are to be encouraged. 

The biological importance of electromechanical coupling 
was discussed in connection with the function of the hair cells 
of the outer ear ( |57ll58] ). In a recent paper, Brownell and col- 
laborators showed that tethers pulled from hair cells contract 
upon application of a voltage difference across the cell body 
(|59l). The phenomena discussed in the present publication are 
especially important for nerve pulse propagation. A change of 
the membrane state from fluid to gel will alter the capacitance 
by about 50%, see eq. (|9]). If this change occurs within 0.5 
ms (which is the time scale of the rising phase of the nervous 
impulse), it will generate a significant capacitive current of or- 
der 100 /iA/cm^ at a constant voltage of Vm = 100 mV. Ionic 
currents of a similar order of magnitude are central elements 
in the Hodgkin-Huxley model of the nerve pulse (3) (Fig. 18 
therein). Hodgkin and Huxley calculated net ion currents of the 
order of 100-600 /iA/cm^ in the squid axon. In previous pub- 
lications we have suggested that the nervous impulse consists 
of an electromechanical soliton corresponding to a lateral com- 
pression of the neuronal membranes (13, 18 -22 ). In particular, 
we have proposed that the membrane undergoes a change in 
state from fluid to gel and back during the nerve pulse. Thus, 
in the soliton theory one expects capacitive currents of similar 
magnitude as the ionic currents in the Hodgkin-Huxley model. 
If charges cannot dissociate (i.e., no capacitive current), one 
rather expects changes in transmembrane potential (Fig. |5]). 
The consequence of these non-linear effects is an electrome- 
chanical soliton that can travel along membrane cylinder with 



Conclusion 

We have shown here that lipid membranes are very susceptible 
to voltage changes close to phase transitions. This result ex- 
tends previous theoretical and experimental findings of voltage- 
dependent capacitances of artificial membranes far from such 
transitions. We have introduced a capacitive susceptibility that 
displays a pronounced maximum at the transition. Voltage chan- 
ges generate changes in area that result in an electromechanical 
coupling. Since biomembranes are exist naturally in a state 
close to a transition, this effect will play a role in excitatory 
processes of the cell. One important example is nerve pulse 
propagation. 
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